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Abstract 

The contractor REnormalization group (CORE) approximation, a new method for solving 
Hamiltonian lattice systems, is introduced. The approach combines variational and contraction 
techniques with the real-space renormalization group approach and is systematically improvable. 
Since it applies to lattice systems of infinite extent, the method is suitable for studying critical 
phenomena and phase structure; systems with dynamical fermions can also be treated. The method 
is tested using the 1+1-dimensional Ising model. 



1^ 

. Introduction 

o 

O^erturbative methods are inadequate for investigating 
^5nany current problems in high energy and condensed 
Q^natter physics, such as the confinement of quarks 
"j^nd gluons in QCD. Thus, the development of new 
C^onperturbative tools is important. The COntractor 
|jk,Enormalization group (CORE) approximation, a 
C^w nonperturbative method for studying Hamiltonian 
^^ ttice systems, is presented in this talk. The method 
IS a hybrid of contraction, variational, cluster, mean- 
field, and block renormalization-group techniques. It 
is systematically improvable and applies to lattice 
systems of infinite extent, enabling direct study of phase 
structure and critical phenomena. Dynamical fermions 
can be treated without problem. 

We briefly describe the method, then apply two 
variants of the CORE approximation to the 1+1- 
dimensional Ising model. 

2. Description of the Method 

The success of any variational calculation, especially 
one involving an infinite number of degrees of freedom, 
depends crucially on choosing a good trial state. An 
algorithm for building trial states suitable for lattice 
systems is the Hamiltonian real-space renormalization 



group (RSRG) method [1]. In this approach, the lattice 
is partitioned into blocks including a few sites and the 
block-Hamiltonians are diagonalized. The Hilbert space 
is then thinned, by discarding all high-energy states, 
retaining only those states which can be constructed 
from tensor products of some small subset of low- 
lying block eigenstates, and an effective Hamiltonian 
which describes the mixing of the remaining states is 
computed. This thinning process is repeated again and 
again until the effective Hamiltonian takes a fixed form 
which can be diagonalized. 

Unfortunately, simple RSRG truncation procedures 
often have difficulties accurately describing the long- 
wavelength modes on the full lattice because they badly 
underestimate the block-to-block mixings. Past ap- 
proaches to overcoming this problem have concentrated 
on using larger blocks, increasing the number of states 
retained per block, or introducing more sophisticated 
truncation schemes. The CORE approximation is a new 
approach to this problem which emphasizes simplicity 
and versatility; it frees one from the need to develop 
clever truncation schemes and allows the use of mani- 
festly gauge-invariant RSRG schemes when studying lat- 
tice gauge theories. 

The basic idea of the CORE approach is to steer 
the RSRG iteration using contraction techniques. An 
important part of this steering process is reliably 
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approximating the expectation value 
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which tends to the lowest eigenvalue eo of Hamiltonian 
H as t —f oo, assuming the trial state I'i'var) has non- 
vanishing overlap with the ground state of H. Since 
e~*^ cannot be computed exactly, it is replaced in 
the CORE method by an operator T(t) which closely 
approximates e~*^ for t in some range < t < tmax 
and which can be easily evaluated. A procedure for 
constructing such an operator has been described in 
Ref. [2]. In this procedure, one expresses e~*^ as a 
symmetric product of explicitly computable terms; for 
example, if iJ = Hi + H2, where e~*^^ and e~*^^ can 
be evaluated exactly, then 
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where C3(t) is order or higher. To construct T(t), 
one then either replaces 6*^^^*) by the identity operator 
or retains low-order terms in Cz{t), rewriting their 
exponential again as a symmetric product of computable 
terms. A given contractor T{t) can also be improved by 

usmgTj,(t) = [T(Vp)f. 

Having chosen a contractor T{t), a variational best 
estimate for eo can be obtained by minimizing 



{^,^r\T{t)HT{t) 1$ 
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with respect to t and any parameters in I'i'var)- For 
a trial state I'i'var) = Yl'j=i '^jl'f'j) 1 where is 
any set of orthonormal states, minimizing Sxit) with 
respect to the aj parameters is equivalent to solving the 
generalized eigenvalue problem 



det {lT{t)HT{t)} - A[T(t)2]) = 0, 



(4) 



where [. . .] denotes truncation to the subspace spanned 
by the \4>j) states. Hence, we can replace the problem of 
finding the best trial state by that of diagonalizing the 
effective Hamiltonian 



ffeffW = [T(t)2]-i/2 lT{t)HT{t)j lT{ty 



-1/2 



(5) 



Developing this operator in the RSRG iteration instead 
of [iJ] is the key innovation of the CORE approach. 

The effective Hamiltonian cannot be exactly deter- 
mined. The last step in the CORE approach is to apply 
cluster techniques to approximate Hf.f[(t) (see Ref. [3] 
and references cited therein). Essentially, this involves 
evaluating Hf.f[(t) on increasingly-larger, connected sub- 
lattices and using the principle of inclusion-exclusion to 
appropriately combine the results for the full lattice. 

In summary, CORE is an iterative blocking and 
thinning process, developing the low-lying physics in a 
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Figure 1. Fractional error 6eo iii the ground-state energy 
density estimates against A. Results using (dashed curve), 
T^^ (solid), and T2 ^ (diamonds with dotted curve) are shown. 

sequence of effective sub-Hamiltonians H^^\t^) using 
the recursion relation 

Hi^+'\t) = R4t)lT(-\t)Hi^\t:)T(-\t)jR^{t), (6) 

where Rn{t) = [T(")(t)2 ]-i/2_ the contractor T(")(t) 
approximates exp[—tH^^\t^)], and is a best value 
for t selected for each RG iteration in some manner: 
minimizing H^^\t) in a simple product state is one 
possibility. As the recursion proceeds, the effective 
Hamiltonian evolves eventually into a simple form which 
can be easily diagonalized, yielding estimates of the 
ground state energy and the energies of some low-lying 
excited states. 

CORE can also be used to estimate the vacuum 
expectation value of an extensive operator O. Using 
the same RSRG transformations as for H, one first 
computes the sequence of effective operators 0^^(t^). 
Once iJeff has evolved sufficiently such that its ground 
state can be found, the matrix element of Oeff in the 
ground state of iJgff then yields the desired expectation 
value. 

3. The l-|-l-Dimensional Ising Model 

The Ising model in 1 + 1 dimensions is often used as 
a testing ground for new calculational methods. Its 
Hamiltonian is given by 
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[cxCT^j) + SX(7x{j)(7x{j + 1)] , (7) 



where j labels the sites in the infinite chain, cx = 
cos(A7r/2), and sa = sin(A7r/2), for 0<A<1. A second- 
order phase transition occurs in this model at A = 1/2. 
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Figure 2. Mass gap estimates A against A. The diamonds and 
squares indicate CORE estimates obtained using T^^ and T2^, 
respectively. The exact mass gap appears as a sohd curve. 



Figure 3. Magnetization against A. The diamonds indicate 
CORE estimates obtained using T2^, the sohd curve shows the 
exact magnetization, and the dot-dashed curve shows the 
estimates from mean-field theory. 



When A < 1/2, the order parameter {ax{j)) = and 
the ground state is unique. For A > 1/2, spontaneous 
symmetry breaking occurs and the order parameter 
takes values (<7^(j)) = ±(1 - cot2(A7r/2))i/®. 

We tested the CORE approximation in two differ- 
ent applications to the Ising model. In both appli- 
cations, the Hilbert space was thinned to the lowest 
two eigenstates in each block, the cluster expansion of 
Hef[(t) was truncated after three-block clusters, and t 
was fixed by minimizing the expectation value of ffgff 
in a mean-field state. Two-site blocking was used in 
the first application, and blocks containing three sites 
were used in the second application. The contrac- 
tor for the first application was Ti(t) = Sl{t)Si{t), 
with Si{t) = YlaiUj [1 + tanh(c„V2)0a(i)]}, where c„ 
are couplings, a labels the different types of operators 
Oa{j), such as cr^i^j) and crj,(j)(7j,(j+l), and j is a site la- 
bel. The second contractor used was T2{t) = Sl{t)S2{t) 
with S2(t) = exp(— ty/2) exp(— tiJj/2), where iJj con- 
tains all intra-block interactions and V contains all 
inter-block operators (those which cross block bound- 
aries). Note that exp(-tiJi,/2) = H exp(-tiJi,(p)/2) 
and exp(—tV/2) = J^^^ exp(— ty(p)/2), where p labels 
the blocks. Calculations were done using T"(t/n) and 
T^(t/n) for various values of n. 

Fractional errors 5^^ = |(£^o — ^o)/^o\ in the ground- 
state energy estimates Eq from both variants of the 
CORE approach are shown in Fig. 1. Selected estimates 
for the mass gap A and magnetization A4 = |(ca;(i))|, 
for some site j, are compared to the exactly-known 
results in Figs. 2 and 3. Considering that only the first 
three terms in the cluster expansion are included in the 
calculations, the accuracy of the results is striking. The 



CORE approximation reproduces the correct location of 
the critical point with remarkable precision. Including 
more terms in the cluster expansion should significantly 
improve these results. 

4. Conclusion 

We believe that the CORE approximation will prove to 
be a powerful tool for studying nonperturbative systems. 
An exciting feature of the method is that it can be 
used to analyze systems containing dynamical fermions, 
systems which resist treatment by present stochastic 
means. We are presently extending the method for use 
with lattice field theories. 
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